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Abstract. We give a mathematical theory of the wetting phenomenon in the 2D Ising model using 
the formalism of Gibbs states. We treat the grand canonical and canonical ensembles. 



1 Introduction 



We study the wetting phenomenon in the 2D Ising model, starting from basic principles of 
Statistical Mechanics. The results of section |3| are based on ]nj, JTTJ] and [|T^] and those of 
section |5] follow from recent results on the large deviations of the magnetization |26|. We 
shall in general refer to these papers for proofs. Our purpose is to give a global view of the 
mathematical results, which are now fairly complete. Since the results about large deviations 
are valid in the 2D case only we restrict the whole discussion to this case. 
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Let us suppose that we have a binary mixture and that the physical parameters are 
chosen so that we have coexistence of the two phases, called + phase and — phase. The 
system is inside a box; the horizontal bottom wall w of the box adsorbs preferentially the 
— phase. If we prepare the system in the + phase we may observe the formation of a thin 
film of the — phase between the wall w and the + phase in the bulk, so that the + phase 
cannot be in contact with this wall. This is the phenomenon of complete wetting of the wall. 
Notice that the total amount of the — phase is not fixed a priori. The appropriate statistical 
ensemble to describe this situation is a grand canonical ensemble. The surface tension f (6) 
is the surface free energy due to an interface between the two coexisting phases, making an 
angle 9 with an horizontal reference line. The contribution to the surface free energy due to 
the wall w when the bulk phase is the + phase is r + . In the case of complete wetting of w 
we expect that r + can be decomposed into r~ + f(0), where r~ is the surface free energy 
due to the wall w in presence of the — phase and f (0) is the surface tension of an horizontal 
interface between the film of the — phase and the bulk + phase. On the other hand, if the 
wall adsorbs preferentially the + phase, then we expect that 

t + = t--t(0). (1.1) 

Indeed, if we impose the — phase in the bulk, then we create an interface between the film 
of the + phase near the wall and the — bulk phase. When complete wetting does not occur 
a stability argument shows that one expects a strict inequality 

\t + -t~\ <f(0). (1.2) 

In such a case the bulk phase is in contact with the wall. In other words the state of the 
system near the wall depends on the nature of the bulk phase. Therefore, depending on the 
property of the wall, we have 

\t + -t~\ <f(0), (1.3) 

with equality if and only if complete wetting holds. This fundamental relation has been 
derived in a thermodynamical setting by Cahn ||. In |K|, |~1| and jjTJ] this situation is 
analysed in the Ising model and the fundamental inequality ( |1.3|) is derived directly from the 
microscopic hamiltonian within the standard setting of Statistical Mechanics. The criterion 
for complete wetting, 

\t + -t~\ =f(0), (1.4) 
is interpreted in terms of unicity of the (surface) Gibbs state. 

It is possible to imagine another situation, which is often used to present the phenomenon 
of wetting of a wall. We have again the same system prepared in the + phase, but we put 
a macroscopic droplet of the — phase inside the box. The droplet is attached to the wall 
w if the wall adsorbs preferentially the — phase. The shape of the droplet depends on the 
interactions between the wall and the binary mixture; the shape of the droplet is given by 
the solution of an isoperimetric problem with constraint PS| (Winterbottom's construction). 
Since the total amount of the — phase is macroscopic and fixed, the relevant ensemble is 
a canonical ensemble. The case of complete wetting of the wall w corresponds to the total 
spreading of the droplet of — phase against the wall. The relation between the shape of the 
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droplet and criterion ( |1.3| ) is made through the contact angle 9 between the droplet and the 
wall (see Fig. [7]), 

cos(#) ■ f(9) - sin(e) ■ = r+ - r~ , (1.5) 

dO 

the so-called Herring- Young equation. A derivation of ( |1.5[ ) is given in section [|. Let us con- 
sider the following two extreme (degenerate) situations. When the wall adsorbs preferentially 
the + phase, then 

r + -t~ = -t(0) < 0. (1.6) 



This corresponds in ( |1.5| ) to 9 = ir, which means that the droplet of — phase is not attached 
to the wall. On the other hand, when the wall adsorbs preferentially the — phase, then 
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f(0) > 



;i-7) 



which corresponds in (|1.5| ) to the case 9 = 0. This means that the droplet does not exist 
as such, but spreads out completely against the wall. We show that the above situations 
can be rigorously derived for the 2D Ising model in an appropriate canonical ensemble. The 
analysis consists in deriving sharp estimates of the large deviations of the magnetization p6 |. 
Our work |26|] extends considerably previous works flT7[| , [fUjfl , [[7|, PHI , [|13] and [|14|] on the 



subject, since we can treat the case of an arbitrary boundary magnetic field. 



2 Ising model 

2.1 Gibbs states 

The lattice is 

Z 2 :={i=(t(l),t(2)): t(i)GZ}. (2.1) 

A spin configuration is a function u defined on Z 2 , t i— > tu(t), with uj[t) = ±1. The Ising 
variable at t is 

a{t){uj) := uj{t) . (2.2) 

An edge of Z 2 , (t,f), is a pair of nearest neighbours sites of the lattice Z 2 . We also call 
edge the unit-length segment in H 2 with end-points t,t'. For each edge (t,f) we have a 
coupling constant J((t,t')) > 0. For each finite subset Ac Z 2 the energy is 

Ha:=~ E J((t,t'))a(t)a(t') . (2.3) 
(t,t')nA^0 

Let uj* be given; the Gibbs measure in A with boundary condition u* and inverse temperature 
(3 is the probability measure 

, um >P(, A ■= / ^(A)- 1 exp(-/3if A (cj)) if u(t) = u*(t) for all t ^ A 



1 otherwise. 
The normalization constant is called partition function. 
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Assume that all coupling constants are equal to one. If we choose u* such that uj*(t) = 1 
for all t, then there is a limiting measure 



lim/^. (2.5) 

The same is true for u*, such that uj*(t) = — 1 for all t. The limiting measure is . These 
two Gibbs states are translation invariant and extremal. There is a unique Gibbs measure if 
and only if = fi~^. This happens if and only if (3 < /3 C , where (3 C is the inverse critical 
temperature. The inverse critical temperature is characterized by the property that there is 
a positive spontaneous magnetization, 

m*(J3) ■ = J a t (u)ii + 'P(du) > (2.6) 

if and only if (3 > (3 C . 



2.2 Contours 



Let uo* be a fixed boundary condition. The usual way of describing the configurations of the 
model is to specify the pairs of nearest neighbours sites e = (t, t') such that a(t)a(t') = — 1. 
Equivalently we specify the dual edges e*, that is the edges of the dual lattice (Z 2 )*, 

(Z 2 )* := { t = (f (1), f (2)) : t(i) + 1/2 G Z } (2.7) 

which cross the edges e. We decompose the set formed by all these dual edges into connected 
components. In |26| we further decompose the connected components into a set of paths, 
called contours, using the rule given in Figure |1|; for details see P6 |. As set of edges the 



paths are disjoint two by two. Some paths are closed and are called closed contours; some 
are open and are called open contours. 
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Deformation rule 



Figure 1: Two connected components giving rise to four contours 



3 Grand canonical ensemble, Cahn's criterion 



3.1 Surface Gibbs states 



Since we are interested in boundary effects, we consider the Ising model on the following 
rectangular box A L (r 1 ,r 2 ) = A L . Let r 1 ,r 2 G IN; we set 

A^(r i, r 2 ) := {t e Z 2 : —r\L < t(l) <nL; < t(2) < 2r 2 L} . (3.1) 

The coupling constants of the model are 

jru ._ / h if t — ( s ; — 1) an d t' = (s, 0) with s E Z, ^ 2) 
1 1 otherwise. 

We consider four different boundary conditions 

,,*(+\ — I a a = ±1, if t(2) < 0, , . 

w W-\6 6 = ±1, if t(2) > 0. (3 ' 3) 

The hamiltonian if^ = H\ L with (a, 6) boundary condition can be written 

H L = - E <r{*W)- E fto^*)- E h ^)- ( 3 - 4 ) 



(t,t'>: <eA L : t£A L :t(2)=2r 2 L 

t,t'eA L t(2)=0 ort(l)=±ni 



The term 



- ^ /ia<7(t) (3.5) 

t(2)=0 

describes the interaction of the binary mixture inside with the bottom wall of the box A^, 
which plays the role of the wall w. The wall adsorbs preferentially the a phase. We interpret 
ha as a real-valued boundary magnetic field and we refer to the 6-part of the boundary 
condition, 

E M*), ( 3 - 6 ) 

t£A L :t(2)=2r 2 L 
or t(l)=±r 1 L 

as the boundary condition. Expectation value with respect to the Gibbs measure in A^ is 
denoted by (-) L ((3,h). When h > and 6 = 1, then all contours of a configuration are 
closed. On the other hand, when h < and 6 = 1, then there is exactly one open contour^ 
in each configuration, denoted below by A, with end-points t^L = {—T\L — 1/2, —1/2) and 
t r> L = {r\L — 1/2, —1/2). The boundary condition ( |3.6| ) specifies the type 6 of the phase in 
the bulk of the system in the following sense 0. Let (3 > (3 C and 5 > 0. Given e > there 

Although there is a close connection between the open contour A and the interface, which is created by 
the bulk + phase and the wall which adsorbs preferentially the — phase, one must not identify the open 
contour with the interface. 

2 For positive h this is proven in |0|; the result is already valid if we replace the condition t(2) > L 1 / 2+l5 



by t(2) > d. The same result holds for h > —h* defined in subsection 3.3 using Lemma 7.1 of |2q ]. We need 
condition i(2) > _L 1 / 2 +' 5 only in the case of complete wetting of the wall by the — phase. The proof is based 
on the following simple fact. The surface tension 9(9) is smooth and has a minimum at 8 = 0; If 8 = _L _1 / 2+<5 , 
then L?(8) - Lt{0)+L/2t"{0)L~ 1+2S = L?(0) + 1/2?"(0)L 2S . Lemma 5.5 in [|| implies that the probability 
that the open contour A visits sites t, t(2) > L 1 / 2+<5 , goes to zero faster than L 2 exp{— 0(L 2S )}. 
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exist L and d such that, VL > L , 

\(a(t))t((3,h)-m*(P)\<e, 
for all t = (t(l),£(2)) G satisfying the conditions 

\t(l)-nL\>d , \t(2) - 2r 2 L\ > d , t(2) > L l/2+s . 



(3.7) 



(3.J 



The surface Gibbs states are the limiting measures 

(■> b (/3,/>) := lim <->*09,fc) 

L— >oo 



(3.9) 



or limiting measures defined by choosing different boundary conditions. The existence of the 



limits when b = ±, as well as the following properties are proven in JlO 



1. The two states ( • ) b (/3, h), b = ±, are extremal Gibbs states. 

2. They are invariant under the translations x ^ x + y, y = (y(l),0). 

3. There is a unique surface Gibbs measure if and only if (-) + (/3,h) = (■ fa). 

3.2 Surface tension and surface free energies 

We define the basic thermodynamic quantities which enter in the description of the wetting 
phenomenon. More details on the surface tension are given in [23] and f24| section 6. Surface 



free energies are studied in |10| and [jTTJ; see also p6| . The definitions we use are standard. 



Explanations for them are given in ||23|| . These definitions are not very satisfactory from a 
conceptual point of view (see beginning of the introduction of [0), since they are defined 
without using a precise notion of interface. This weak point is also a strong point, because 
the notion of interface is a delicate notion, which is difficult to analyse. 



3.2.1 Surface tension 



For latter purposes (see section ^) it is more convenient to parametrize surface tension using 
the normal vector to the interface instead of the angle 9. We consider the model with 
coupling constants equal to one on the whole lattice. Let be the square box 



tt L := {t E Z 2 : —L < t(l), t(2) < L) . 



(3.10) 



Let n = (n(l), n(2)) be a unit vector in 1R 2 and C(n) a line perpendicular to n, passing 
through (0,0). We define a boundary condition uj n by setting (see Fig. ||) 



u n {t) :- 



1 if t is above or on the line C(n), 
— 1 if t is below the line C(n). 



(3.11) 
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Figure 2: Definition of the oj n boundary condition 

The corresponding partition function is Z u "(Oi). By definition, the surface tension of 
an interface perpendicular to n is 

r(n) : = lim -— In ^ ^ . (3.12) 

We usually do not write explicitly the /3-dependence of the surface tension. We extend 
the function f to IR 2 as a positively homogeneous function. By GKS-inequalities it follows 
that the function is subadditive, 

r{x + y)<T{x)+r{y). (3.13) 



It can be shown that f is strictly positive (for all x ^ 0) if and only if j3 > (3 C ||16|| . (This is 
also true for dimensions higher than two.) 

Remark: In dimension two the surface tension is related by duality to the decay-rate 
of the two-point function at the dual inverse temperature [Q, P3 |. The decay-rate of the 



two-point function at the dual inverse temperature can be computed explicitly |19| . The 
surface tension has the following symmetry properties, 

?(n) = r(—n) and f (n) = f(m) if nlm. (3-14) 

It is a smooth function for any (3 < oo. 



3.2.2 Surface free energies 

The coupling constants are given by ( ^.2| ). Free energies are defined up to an arbitrary 
constant. We need only to define r _ — r + , the difference of the contributions of the wall w to 
the free energy, when the bulk phase is the — phase respectively the + phase. The definition 
is similar to the definition of the surface tension. Let be the box (|3.1|). By definition, 



r"-T+:= lim -J- In . (3.15) 

l^oo 2r x L Z+(A L ) K ' 
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We usually do not write explicitly the (/3, /^-dependence of r — r + . The existence of the 
limit is proven in |L1|, as well as the following results. 



1. For any fi, \t~ - t+| < f((l, 0)). Hence r~ - r+ = for all (3 < (3 C . 

2. We have the symmetry 

r-(h) - r + (h) = -(r-(-h) - T+(-h)) . (3.16) 

3. Let f3 > (3 C and n w := (0,-1). For positive h, T~(h) — T + (h) is a positive concave 
function; for all h > 1 

T-(h)-r + (h) = r(n w ). (3.17) 
These results are not restricted to dimension two. 



Remark: In dimension two the quantity r — r + can be computed P2" |. For h positive it 



is equal to the decay-rate of the boundary two-point function at the dual coupling constants. 



3.3 Cahn's criterion and phase diagram 

In all the section the inverse critical temperature f3 > (3 C is fixed. We study the model in 
the grand canonical ensemble; the coupling constants are given by (|3.2j ) and the states are 
the surface Gibbs states ( ■ ) h of subsection |3.1| . For definiteness we choose the + boundary 
condition, so that we have the + phase as bulk phase. 



From the properties of r — r + , which are stated in subsection |3.2.2j , we have the existence 
of a positive h*{(3) such that (see Fig. ||) 

h*(/3) := inf [h > : t~{(3, h) - r + (/3, h) = ?(n w )(f3)} . (3.18) 

The thermodynamical criterion of Cahn states that there is complete wetting of the wall if 
and only if \h\ > h*{(5). 

Remark: The value of h*(/3) is known in dimension two. It was computed by Abraham 
[0|; h*(j3) is the solution of the equation 

exp{2/?} {cosh 2/5- cosh 2(3h* (13)} = sinh2/?. (3.19) 



In fll2f , see also [p5fl , it is shown that the value of h*([3) can also be obtained from the work 



of McCoy and Wu P . 



To make the connection with the surface Gibbs states we use the identity proven in [] I 
T-(P, h) - r+03, h) = £ ((a(0)> + (/3, h') - (a(0))-((3, h')) dh' . (3.20) 



T — T 




h* 1 h 
Figure 3: r~ — r + as a function of h. 



Therefore 



Tin,, 



T 



-(fj,h*)-r + (f3,h*) 

(a(0)) + (P,W-(a(0))-{J3,h!j)dti. 



For all h > h* we must have 



(a(0)) + ((3,h) = (a(0))-((3,h) 



(3.21; 



(3.22) 



since for positive hi 

(<r(0)) + ([3,ti)>{a(0))-({3,ti). (3.23) 

On the other hand, if < h < h*, the concavity of T~(/3,h) — T + (i3,h) as function of h 
implies that 

(a(0)) + (P,h)>(a(0))-((3,h). (3.24) 

These results, together with the properties of the surface Gibbs states mentionned in sub- 
section |3.1| , imply the following interpretation of the Cahn's criterion for complete wetting 
(see also the corresponding phase diagram of Fig. |j). 



1. There is complete wetting of the wall w if and only if there is a unique surface Gibbs 
state. 



2. There is a unique surface Gibbs state if and only if 

\ T -(J3,h)-T + (j3,h)\ =?(n w )(j3) >0. 



(3.25) 



3. The value h* of the boundary magnetic field where the wetting transition takes place 
can be defined as 

h*(P) := inf {( • >+(A h) = (- HA h)} . (3.26) 
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T 



Non-uniqueness of 
surface Gibbs state 



-1 1 h 

Figure 4: Phase diagram. The region of non-uniqueness of the surface Gibbs state is shaded. In 
the other region, there is a single surface Gibbs state. 



Remarks: 1. These results are proven in flTT| ; they are not restricted to dimension two. 

2. In dimension two we have precise information about the behaviour of the open contour 
A entering in the description of the configurations when we have — boundary condition and 
< h < h*. The contours sticks to the wall. There exists a constant K, such that the 
probability that the open contour A does not visit the segment I = {t : t{2) — 0, x\ < t(l) < 
£2} is smaller than 

K exp{— k ■ \xi — x 2 \} , (3.27) 

with 

« = t(0 - (t- - r+) . (3.28) 

This result is a combination of Lemma 7.1 in p6[ and the remark following the proof of 
that lemma. By symmetry the same result applies to the case where we have + boundary 
condition and > h > —h*. 

3. Weaker results about the behaviour of A, but not restricted to dimension two, are 



contained in 11 



4 Variational problem 

We consider here the variational problem giving the shape of the macroscopic droplet when 
the surface tension and the surface free energies of the wall are known. In section [5] we show 
how this variational problem arises when the analysis starts from the hamiltonian of the 
model. 

The variational problem, which gives the shape of the macroscopic droplet in presence of 
a wall, is a generalization of the classical isoperimetric problem. In the physics literature the 
solution of the problem is known as Winterbottom's construction J29j. Wulff's construction 



I corresponds to the special case when the wall has no effect on the droplet. Dinghas 
gave a geometrical proof of Wulff's construction, which has been extended by Taylor 
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(see her review for original references). Wulff's solution and Bonnesen's inequalities, 
which describe a (strong) stability property of the solution, are discussed in details in [[?]. 
A completely different proof, valid in the 2D case only, is given in |J[]. The many facets of 
statistical mechanics of equilibrium shapes are reviewed in [27| and [ |3Tf . 



4.1 Geometry of the boundary of a convex body 

Let C C 1R 2 be a compact convex body, that is, a compact convex set with a non-empty 
interior. We denote by dC the boundary of C and by (x\y) the Euclidean scalar product of 
x,y G 1R 2 . The support function of C, Tc, is defined on 1R 2 by 

r c (y) := sup(x\y) . (4.1) 



It is immediate from the definition ( |3~T ) that tq is positively homogeneous and subadditive, 



T C {vx + y%) < r c {yi) + T C (y 2 ) ■ (4.2) 

Since C is compact it is also clear that for any y there exists x y G C, in fact x y G dC, with 

r c (y) = (x y \y) . (4.3) 

Let n G 1R 2 be of norm ||n|| = 1. Define A(n) to be the hyperplane 

A(n) := {x G 1R 2 : (x\n) = T C {n)} . (4.4) 

Then A(n) is a support plane for C at x n , that is x n G A(n) and 

C C {x G 1R 2 : (x\n) < r c (n)} . (4.5) 

Conversely, if x G dC, then there exists a support plane for C at x. This is a consequence of 
the separation result: if O is an open convex set and L an affine subset, such that OC\L = 0, 
then there exists a hyperplane if with the properties || 

Lcif and OCiH = . (4.6) 

For n G 1R 2 , ||n|| = 1, Tc(n) gives the (signed) distance of the hyperplane A(n) to the origin. 
The distance is positive if and only if G {x G H 2 : (x\n) < Tc(n)}. 

Theorem 4.1 Let C be a compact convex body in 1R 2 . Let tq be its support function. Then 

C = {x G 1R 2 : (x\n) < r c (n) Vn, ||n|| = 1} (4.7) 
= p| {x G IR 2 : (x|n) < r c (n)} . 

n: ||'/(|| = 1 



3 These few references are far from complete. Good reviews about the isoperimetric problem and related 
topics are [EOl and [Elf. 
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Proof. It is clear that 

C C p| {x G IR 2 : (x\n) < r c (n)} . (4.8) 

n: ||n||=l 

Suppose that y G" C. We can separate strictly a closed convex set B and a compact convex 
set K by a hyperplane, when they are disjoint ||. Therefore there exists a hyperplane 

if = {x G IR 2 : (x|m) = 5} , (4.9) 

||m|| = 1, such that for all x G C we have (x\m) < 5 and at the same time (y\m) > 5. 
Therefore 

sup(x|m) = Tcijn) < 5 , (4-10) 

and consequently 

y£ p| {x EM 2 : (x\n) <r c {n)} . (4.11) 

n: ||n||=l 

□ 



We recall two definitions (see Fig. |5|). A point y G <9C is a regular point if there is 
a single support plane containing y. This is equivalent to say that the intersection of all 
support planes containing y is a 1-dimensional afline set. At a regular point x there is a 
well-defined (unit) normal vector to C, defined as the vector n such that x G A(ri). A 
support plane A is a regular support plane if A n C is 0-dimensional. 



Theorem 4.2 Lei C be a compact convex body. Let n G H 2 , ||n|| = I. Then A{n) D C 
coincides with the subdifferential drc{n) of tc at n, that is, 

A(n) DC = {x G H 2 : r c (z + n) > r c (n) + G 1R 2 } . (4.12) 

Proof. Suppose that x G drc(n), 

t c (z + n) > r c (n) + (z\x) Vz G H 2 . (4.13) 

Since tc is subadditive, 

t C (z) + r c (n) > r c (n) + (z\x) , (4.14) 

we have 

tc(z)>(z\x) Vz; (4.15) 



hence x G C by Theorem |4.1| . On the other hand, 

(z + n\x) — r c (z + n) < (n\x) — r c {n) ; (4.16) 

therefore 



(x) = sup f (y\x) - T C (y) } = (n\x) - r c {n) . (4.17) 
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Figure 5: Convex body C with support planes A(m), A{n2)- £2 is a regular point and A{n2) is a 
regular support plane. x\ is a regular point and A(n\) is not regular. The vectors y\ and 7/2 are in 
the subdifferential of r at n\. dist(0,Pi) = r(ni); dist(0,i-2) = T{n-i). 

Since x G C, we have for all y 

(y|x)-r c (y) <0; (4.18) 

hence 

sup {(y\x)-T C (y)} = 0, (4.19) 

that is x G A(n). 

Conversely, suppose that x G A{n) n C. Since x G we have 

(x|n) = T C {n) . (4.20) 

Since x G C we have for all z 

r c (^ + n) > (z + n\x) . (4.21) 

Therefore 

r c (^ + n) > (n\x) + (z\n) = r c {n) + (z\n) . (4.22) 

□ 

As a consequence of Theorem [O the support plane A(n) is regular if and only if tq is 



differentiate at n, \\n\\ = 1. On the other hand the regular points of dC are characterized 
by the points of dC where the normal is well-defined (there is a unique tangent plane at 
those points). We can partition the boundary dC of C into three sets. 



1. x is regular and the support plane A(n), n the unit normal vector to C at x, is regular. 
The support function t c is different iable at n. 
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2. x is regular and the support plane A(n), n the unit normal vector to C at x, is not 
regular. The support function tq is not different iable at n. The facet of (unit) normal 
n is given by the subdifferential of the support function tq at n, 

A{n) nC = dr c . (4.23) 

3. x is not regular. Such a point corresponds to a corner of C, where the normal is not 
well-defined. There are at least two different support planes. 



Remarks: 1. The dual function of tq is 

r* (x) := sup { (x\y) - r c (y) } . (4.24) 
yen 2 L J 

In our case the dual function r£ is the indicator function of C: 

(See proof of Theorem |4.2| .) The dual function of Tq is tc, 

t**(x) := sup { (x\y) - r* c (y) } = r c (x) . (4.26) 
yen 2 L J 

2. There is another interesting function associated to C, the polar of C, Tq, 

t°c(x):= sup £W (4.27) 

ne]R 2 ||n||=l Tclnj 

The convex body C can be expressed as 

C = {x G 1R 2 : rc(x) < l} , (4.28) 

and its boundary as 

DC = [x G H 2 : ro(ar) = l} . (4.29) 
A point of DC is regular if and only if Tq is differentiable at x. 



4.2 Isoperimetric inequality 

Let C be a compact convex body in H 2 and tq its support function. Let dV be a rectifiable 
curve in 1R 2 , which is the boundary of an open set V. The curve dV is oriented in the usual 
way; n(x) denotes the exterior unit normal vector at a point x G dV. We define a functional 

F c (dV) := / r c (n(s))ds. (4.30) 

JdV 
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If 7 : [0, t] — > IR 2 is a parametrization of the boundary dV, then, using the homogeneity 
property of tc, we can write the functional as 

T c {dV) := jT* r c ((7(2)(tz),-7(l)(w)))d«, (4.31) 



where 7(5) is the derivative of 7 with respect to s. The functional is always positive. Indeed, 
by a suitable translation a we can suppose that the translated set C = C + a has zero as 
interior point, and consequently its support function tc> is strictly positive at x 7^ 0. Since 

T C '(x)=T C (x) + (x\a), (4.32) 

and 

/ (n(s)\a)ds = 0, (4.33) 
the value of the functional does not depend on a. 

Theorem 4.3 (Generalized isoperimetric inequality) Let C be a compact convex body 
in 1R 2 . Let V be an open set in IR 2 such that its boundary is a rectifiablef\ curve dV . The 
Lebesgue measure of C and V is denoted by \C\ and \ V\. Then 

Fc(dV) >2\C\ 1/2 \V\ 1/2 . (4.34) 

Equality holds in ( 4-3$) if and only ifV equals, up to dilation and translation, the set C . 



Remarks: 1. The set V may have several connected components. 

2. The main ingredient of Dinghas' proof is to express the functional Tc in a geomet- 
rical manner as 

\V + eC\ - \V\ 

F c (dV) = lim 1 - 1 ^ , (4.35) 

slO E 

where V + eC denotes the set 

IJ (x + eC) , (4.36) 

and eC = {ex : x G C}. The isoperimetric inequality follows by applying Brunn-Minkowski 
inequality to V + eC, 

\V + eC\ 1/2 > \V\ 1/2 + \eC\ 1/2 . (4.37) 

3. When C is the unit ball we recover the classical isoperimetric inequality. 

4. The minimum of the functional for open sets V with \V\ = \C\ can easily be computed, 

min FMdV) = 2 • |C| . (4.38) 

V:|V|=|q 



The theorem is valid under less restrictive assumptions; see e.g. 



15 



In addition to the isoperimetric inequality, the stability of the minimum can be controlled 
in a rather strong sense by the (generalized) Bonnesen's inequalities. Let 



and 
Then 



r(dV) = sup jr : r ■ C + x C V for some x G 1R 2 | 
R(dV) = inf [R : i? • C + x D V for some x G H 2 } 



T c {dV)- T c {dVf-4\C\-\V\ 



1/2 



2|C| 



< r{dV) < R{dV) 



(4.39) 
(4.40) 

(4.41) 



< 



Fc{dV) + 


>o(^) 2 - 4|(7| • 


1/2 


2|C 





4.3 Variational problem 

Let ri,r 2 G IN be given and define the rectangle Q, 

Q--{ x = ( x (l),x(2)) G H 2 : -r x < x{l) < n ; < x(2) < 2r 2 } . (4.42) 

Let (3 > (3 C and f : 1R 2 — > 1R be the surface tension. Let S 1 = {n £ 1R 2 : ||n|| = 1}; we define 
on Q x S 1 a function t, 

/ v f r _ — r + if x(2) = and n = n w , , . 

1 f(n) otherwise. 

We define a functional as above, by setting 

W t (dV) := / t(s;n(s))ds. (4.44) 

Notice that r _ — r + is negative when /i is negative. However, as long as r~ — r + > -f(n w ), 
then the functional is positive. This is an easy consequence of the elementary monotonicity 



principle stated in subsection [4.3.1 



Let us give the interpretation of W t (7) in our setting. We consider a macroscopic droplet 
of the — phase immersed in the + phase. The boundary of the droplet is dV. Suppose first 
that x(2) ^ and ||n|| = 1. Then t(x; n) is interpreted as the surface tension of an interface 
perpendicular to n, passing through x. The case x{2) =0 and n = n w corresponds to the 
situation where the macroscopic droplet is in contact with the wall. In that case 

t- - t + (4.45) 

is the change in the free energy due to the presence of the droplet on the wall (the bulk 
phase is the + phase). 
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Variational problem VP: Suppose that the bulk phase in Q is the + phase and that 
the — phase occupies a set V such that \V\ is a fraction a of the volume of Q, \V\ = a\Q\. 
V is not necessarily connected, but we assume that it is open and that its boundary dV is a 
rectihable curve. Find the optimal set C , \C\ = ct\Q\, such that C minimizes the functional 
V t (dV),\V\ = a\Q\. 



Remarks: 1. In |26j we consider the same variational problem. Its formulation is 
slightly different, because we have introduced as fundamental quantities the dual quantities, 
the decay-rates of the two-point functions. 

2. There is a simple way of showing that the solution C is a convex body. Let us consider 
the case where t(x,n) = t(ti). Since f is convex and positively homogeneous, we have for 
any parametrized curve 7 : [a, b] — > 1R 2 , by Jensen's inequality 



f((7(2)(s),-7(l)(s)))ds> (4.46) 



"(( 



b — a 

7(2)(6)- 7 (2)(o) 7(1)(&) -7(l)(a 



b — a b — a 

Therefore one decreases the value of the functional each time we replace some part of the 
curve between two points by the straight segment between these two points. Since the box 
Q is convex, for every open set V C Q its convex enveloppe conv V C Q, and |conv V\ > \V\. 
On the other hand 

W t (OV) > W t (<9(conv V)) . (4.47) 



4.3.1 Solution of the variational problem 

We give the explicit solution in some special cases. We refer to [TI| for other cases. We 
consider cases where only the effect of the bottom horizontal wall is important. This amounts 
to consider cases where the amount of — phase is not too large. (Theorem |5.1| is valid without 
such restrictions.) 

The function f can be interpreted as the support function of a compact convex set , which 
we denote by C9, 

CV:={xelR 2 : (x\y)<T(y)Vy}. (4.48) 

Indeed, the dual function f* is the indicator function of the set C9, so that the support 
function of C~, which the dual of the indicator function is r ** = f, because f is convex. The 
compact convex body C~ is called Wulff shape in the physics literature. 

Case 1: r~ — t + = f(n w ). 

In that case t is independent of the variable x, so that 

t(x,n) = f(n) . (4.49) 
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Figure 6: Winterbottom construction. 



Let us fix the volume of the set V occupied by the — phase tobea|Q|,0<a<l. Ignoring 
for a moment the constraint that V C Q, the solution of the variational problem is given by 



Theorem fj~3| . It is a Wulff shape of volume a\Q\, A • C~, with 



A=g. (4.50) 

Whenever a translate of A • can be put inside Q, this is a solution of the variational 
problem. When there is no translate of A • C9 which can be put inside Q, then the constraint 
that V C Q modifies the shape of the optimal set V (see []15||). 

Case 2: |t~ — t + | < f(n w ). 

Let us first ignore the constraint that V G Q. Notice that the problem without the 
constraint V C Q is scale-invariant. Then the solution is given by a Winterbottom shape of 
volume a\Q\. The Winterbottom shape is by definition the convex body (see Fig. ^) 

C t :=Cf]{x: (x\n w ) < r~ - r + } . (4.51) 



To prove the optimality of the Winterbottom shape we use Theorem |4.3| and the following 
monotonicity principle Suppose that we can find a convex body C, such that the 
following two conditions ( |4.52| ) and (|4.53|) are verified. If we replace in the definition of the 
functional W t the function t by the support function Tc, then 



and for V — C, 
Then we have 
and 



Mt(dV)>V 7V (dV), 

W t (0C) =W TC (<9C). 
MdV) > 2\C\ 1/2 \V\ 1/2 



V t (dV)- U t (dV) 2 - 4\C\ ■ \v\ 



1/2 



2\C\ 



< r(dV) < R(dV) 



(4.52) 

(4.53) 
(4.54) 

(4.55) 
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< 



M9V) + 


MdV) 2 


-4\C\- 


v\ 


1/2 


2\C 





The proof of Q4.54| ) is an immediate consequence of Theorem |4.3| applied to the convex set 
C . The proof of (|4.55| ) is an immediate consequence of the monotonicity of the real function 

s^s-(s 2 - D) 1/2 (4.56) 

for \s\ > \f~D. In our case C = C t . It is immediate that condition ( |4.53|) is verified. To 
show ( 4.52 ) we notice that t* is the indicator function of C t . Therefore its support function 
satisfies 

t(n) > t**(n) = T C {n) . (4.57) 



Any translate of the Winterbottom shape, which is contained in Q, is a solution of the 
variational problem. When there is no translate which is inside Q, then the constraint that 
V C Q effectively modifies the shape of the optimal set V. We want to discuss one simple 
situation of that kind (see Fig. || and |9] for illustrations). Suppose that a is chosen so that 
a translate of the Winterbottom shape of volume a\Q\ exists inside Q whenever h > 0. Let 
h(a) be the smallest value of h such that a translate of the Winterbottom shape of volume 
a\Q\ exists inside Q. Notice that necessarily h(a) > —h*, and that the value of h(a) is a 
function of the box Q. Then for any h < h(a) the solution of the problem is the same as the 
solution of the problem for h = h(a): the box prevents the macroscopic droplet to spread 
out. For h < h(a) the optimal shape of the droplet of volume a\Q\ is not the corresponding 
Winterbottom shape. 




Figure 7: Angles 9 and ip 



Remark: The Young-Herring relation, giving the contact angle 6 of the macroscopic 
droplet with the wall in terms of the surface tension and surface free energies, can be easily 
derived from Theorem [4.2| . Let us parametrize the unit vectors in 1R 2 by an angle cp as in 
figure [7[ Let n be the normal to the interface defining the contact angle 9. Theorem ^L2 
asserts that 

gradr(n)(2) = ~(r" - r+) . (4.58) 
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This is exactly Young-Herring relation. Indeed, in polar coordinates the surface tension is 

9 (r, <p) :— r ■ r{n{<p)) , (4.59) 
so that, using x(2) = rsinyj and <p + 9 — ir/2, we get (f(l,<p) = f((p)) 



gradf(ra)(2) = sin (p? (<p) + cosy? —f{(p) 

dip 

= cos 6 t(6) - s in0-^r(0) 
do 



(4.60) 



a 


b 


c 


d 


e 


f 









Figure 8: a: h > h*; b: < h < h*; c: h = 0; d, e, f: a sequence of droplets for decreasing values of 
— h* < h < 0. The droplet spreads until it begins to touch the vertical sides of the box (e). Further 
reduction of the magnetic field does not modify the shape of the droplet, but makes it unstable in 
the sense that the removal of the vertical walls would result in a spreading of the droplet, f: the 
dashed line shows only a part of the droplet which would be obtained by removing the walls when 
> h > -h*. 



5 Canonical ensemble and macroscopic droplet 



To study the macroscopic droplet of fixed volume we introduce a canonical ensemble. There 
is some freedom in doing this, in the sense that we can fix the total magnetization up to 
fluctuations, which are negligible when measured at the scale of the volume. In subsection 



O we define the canonical states and in subsection [5.2| we consider the limit of the lattice 
spacing going to zero. We show that in this limit there is a droplet of the — phase immersed 
in the + phase and that the shape of the droplet is given by the variational problem of 
subsection [4.3| . This completes the mathematical theory of wetting in the 2D Ising model in 
terms of the Gibbs states. We always consider the case of + boundary condition. 
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Figure 9: Sequence of big droplets in a tube for decreasing value of the magnetic field. a,b: The 
upper part of the droplet has the Wulff shape, while the lower part has the winterbottom shape; 
they are joined by a rectangle such that the total volume is conserved; c: The droplet completely 
wets the lower wall, the droplet is build up from a half Wulff shape and a rectangle. This situation 
holds as long as — h* < h < 0. 

5.1 Canonical states 

We define the canonical states at finite volume. Let (3 > j3 c and — m* < m < m*. Let 
c = 1/4 — 6 > 0, with 5 > 0. We introduce the event 



A{m- c) :={<*>: | £ u(t) - m|A L | | < \K L \ ■ L~ c } 
teA L 



(5.1) 



The canonical state in A, with + boundary condition and parameter m, is the condi- 
tional state 

(•\m)i(p,h):={>\A{ m] c))i(p,h). (5.2) 

To understand the canonical state (|5.2|) we must control the large deviations of the magneti- 
zation in the state (•)J I - Although for typical set of configurations with respect to (•)"£ the 
length of the contours is small f], because we impose here a specific magnetization m ^ m* 
there is always at least one large contour in each typical set of configurations of the canonical 
state. It is therefore natural to distinguish between large and small contours. Let 5(0; [L s ]) 
be the square box 

5(0; [L 5 ]) := {teJR 2 : \t(j)\ < [L 5 ] , j = 1,2} . (5.3) 

We say that a contour 7 is small if there exits a translate of the box 5(0; [L 5 ]) which 
contains 7. Otherwise the contour is large. We sum over small contours and the large 

5 For any (3 > (3 C there exits K{0) so that limL^oo ( {3 a contour of length > K\nL] )J = . (Lemma 
5.6 in Hi.) 
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contours are treated by a coarse-graining method. Theorems 11.1 and 11.2 in [T2(J give a 
detailled description of a set of typical configurations in terms of large contours. One result 
of this analysis is the exact computation of the large deviations of the magnetization for the 
Gibbs state ( • h), together with a control of the speed of convergence. 

Theorem 5.1 Assume that (3 > f3 c , h e 1R ; —m*{(3) < m < m*(f3) and c := 1/4 - 8, 8 > 0. 
Let W*(m) be defined by 

s <rff * — 777 

r t (m) := inf {v t (C) : C C Q , volC = -^-—\Q\). (5.4) 
Then for any rj < 5 and L large enough 

| i In P h L [A{m- c)} + w;(m) | < 0(L^ 5 ) ; (5.5) 
the probability is computed with the measure ( ■ )J(/3, h). 



5.2 Macroscopic droplet 

In this last subsection we consider the limit of the lattice spacing going to zero. We suppose 
that (5 > 0c, h 6 M and we choose the + boundary condition. The probability measure in 
this section is always the canonical Gibbs state ( • |m)£(/3, h) defined in subsection |54] with 
— m* < m < m* and c = 1/4 — 5 > fixed. We do the analysis in the box A^ri,^) and 
at the end we scale everything by 1/L, i.e. we take the limit of the lattice spacing going to 
zero. 

Let C C Z 2 ; the empirical magnetization in C is 

mc(w):=|4E *(*)M- M) 
M tec 

Let < a < 1; we introduce a grid C(a) in made of cells which are translates of the 
square box 

5(0; [L a ]) = {t E H 2 : \t(j)\ < [L a ] , j = 1,2} . (5.7) 

The value of a is close to 1. In most of the cells the empirical magnetization is close to 
m* or — m* with high probability. For each cell of the grid C(a) we compute the empirical 
magnetization mc(w). Then we scale all lengths by 1/L, so that after scaling the box is 
the rectangle Q. For each u we define a magnetization profile Pl(%', on Q, 

Pl{x;u) := mc{uj) if Lx G C (5.8) 

where Lx is the point x & Q scaled by L and C a cell of the grid C(a). 
The set of macroscopic droplets at equilibrium is 

777 * — 777 

V(m) := {V C Q : \V\ = ^— -\Q\ , W t (<9V) = W t *(m) } . (5.9) 
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For each V G T>(m) we have a magnetization profile, 

/ \ ( m* if x G Q\V, /riM 
Pv{x) :=\ , c 5 - 10 
L — m x G V. 

Let / be a real-valued function on Q; we set 

d x {f,V{m)):= inf / do; | / (x) - p v (*) i • (5.11) 
The main theorem (Theorem 12.2 in p6| ) is 

Theorem 5.2 Let (3 > (3 C , h G 1R ; -m* < m < m* , c = 1/4 - 5 > 0. Let ( • | m )J(/3, /i) 5e 
t/ie canonical Gibbs state with + boundary condition. Then there exists a positive function 
s(L) such that limi^ 00 e(L) = and for L large enough 

Prob[ { d 1 (p L ( ■ ; w), P(m)) < e(L) }] > 1 - exp{-0(L K )} . (5.12) 

This theorem gives a complete description of the wetting phenomenon in the canonical 
state, making the connection between a microscopic approach with a conditional state and 
the macroscopic variational problem giving the shape of a macroscopic droplet at equilib- 
rium. The two approaches with grand canonical and canonical ensembles give of course the 
same information about the occurence of the wetting transition. However, for this surface 
phenomenon they are not equivalent, but complementary. 
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